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Context. T Tau stars display different X-ray properties depending on whether they are accreting (classical T Tau stars; CTTS) or not (weak-line 
T Tau stars; WTTS). X-ray properties may provide insight into the accretion process between disk and stellar surface. 

Aims. We use data from the XMM-Newton Extended Survey of the Taurus Molecular Cloud (XEST) to study differences in X-ray properties 
between CTTS and WTTS. 

Methods. XEST data are used to perform correlation and regression analysis between X-ray parameters and stellar properties. 
Results. We confirm the existence of a X-ray luminosity (Lx) vs. mass (M) relation, Lx <x M L69±0 U , but this relation is a consequence of 
X-ray saturation and a mass vs. bolometric luminosity (L.) relation for the TTS with an average age of 2.4 Myr. X-ray saturation indicates 
L x = constL,, although the constant is different for the two subsamples: const = lO" 3 - 73 * 005 for CTTS and const = lO" 3 - 39 * 006 for WTTS. 
Given a similar L, distribution of both samples, the X-ray luminosity function also reflects a real X-ray deficiency in CTTS, by a factor of w 2 
compared to WTTS. The average electron temperatures L av are correlated with L x in WTTS but not in CTTS; CTTS sources are on average 
hotter than WTTS sources. At best marginal dependencies are found between X-ray properties and mass accretion rates or age. 
Conclusions. The most fundamental properties are the two saturation laws, indicating suppressed Lx for CTTS. We speculate that some of 
the accreting material in CTTS is cooling active regions to temperatures that may not significantly emit in the X-ray band, and if they do, 
high-resolution spectroscopy may be required to identify lines formed in such plasma, while CCD cameras do not detect these components. 
The similarity of the Lx vs. T m dependencies in WTTS and main-sequence stars as well as their similar X-ray saturation laws suggests similar 
physical processes for the hot plasma, i.e., heating and radiation of a magnetic corona. 

Key words. Stars: coronae - Stars: formation - Stars: pre-main sequence - X-rays: stars 



1. Introduction 

Optically revealed low mass pre-main-sequence stars define the 
class of T Tauri Stars (TTS). TTS are divided into two families, 
the Classical T Tauri Stars (CTTS) and the Weak-line T Tauri 
Stars (WTTS). CTTS display strong Ha lines, a sign that the 
stars are accreting material from the circumstellar disk, while 
in WTTS the Ha line fluxes are suppressed, a sign that accre- 
tion has ceased. Based on infrared observations, Young Stellar 
Objects (YSO) have instead been ordered in classes accord- 
ing to their infrared (IR) excess. Following this classification, 
deeply embedded stars at the start of their accretion phase are 
"Class 0" objects, more evolved protostars still embedded in 
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their envelope are "Class I" objects, stars with a circumstel- 
lar disk that show IR excess are "Class II" objects, and stars 
with no IR excess are "Class III" objects. While the Ha clas- 
sification is based on accretion, the IR excess is a measure of 
circumstellar material. The Class II objects are dominated by 
CTTS, while the Class III stars are dominated by WTTS. 

Both types of TTS have been found to be strong X-ray 
emitters. First X-ray detections of ind ividual TTS were made 
with the Einstein observatory (e.g., IFeigelson & DeCampli , 
19811) and revealed very strong X-ray activity, exceeding 
the solar level by several orders of magnitude. Many star- 
forming regions have s ubsequently been obse r ved with the 



ROSAT satellite (e.g. IFeigelson et all Il993t iGagneetal 



19951: NeuhauseretaU 119951: |Stelzer_& Neuhauser, 2001), 



largely increasing the number of X-ray detected TTS. Studies 
based on Ha emission may in fact fail to detect part of the 
WTTS population, which can easily be identified in X-rays. 
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The origin of the strong X-ray activity in TTS is not en- 
tirely clear. The observed emission in the soft X-ray band above 
1 keV is consistent with emission from a scaled-up version 
of the solar corona. In main-sequence stars X-ray activity is 
mainly determined by the stellar rotation rate. The activity - 
rotation relation is given byLx/L* oc P; 2 - 6 dGiidel et allll997l) . 



where Lx is the X-ray luminosity, L* is the stellar photospheric 
bolometric luminosity, and P rot is the rotation period of the 
star. This is consistent with the dynamo mechanism that is 
present in our Sun, wher e the magnetic fields are generated 
through an a-Q, dynamo dParkei , fl955l) . At rotation periods 
shorter than 2-3 day s for G-K stars, the X-ray activity saturates 
at log(L x /i.) ~ -3 dVilhu&Rucinskilll983h . 



As for pre-main sequence stars , early surveys of the 



Taurus Molecular Cloud (TMC; iNeuhauser et all 11995 



Stel zer & Neuhauserl 120011) claimed a rotation-activity rela- 
tion somewhat similar to the relation for main-sequence stars, 
but the recent COUP survey of the Orion Nebula Cluster 



(ONC ) found absence of such a relation dPreibisch et al 



2005), suggesting that all stars are in a saturation regime, 
even for long rotation periods. Young stellar objects, espe- 
cially in their early evolutionary stage, are thought to be 
fully convective, and the generation of magnetic fields through 
the a-Q dynamo should not be possible. This suggests that 
X-rays in low-mass pre-main sequence stars are generated 
through processes different than in the Sun. New models for 
X-ray gene ration through other dynamo concepts ha v e bee n 
developed ( Kiiker & Riidigei . 19991 : Giampapa et al. . 1996). 
Alternatively, in CTTS, X-rays could in prin ciple be produced 
by magnetic star-disk interactions (e.g., iMontmerle et al. , 
2000[ llsobe et all 120031). i n accretion shocks ( e.g.. [Larnzin. 
1999J iKastner et all l2002t Istelzer & Schmittl l2004 or in 



shocks at the base of outflows and jets dGiidel et all 12005 



Kast ner et all 12005). 



The influence of a circumstellar disk, and particularly the 
influence of accretion on X-ray activity is therefore of in- 
terest. Former X-ray studies of star forming regions have 
led to discrepant results. In the Taurus-Auriga complex, 
Stelzer & Neuhauserl d200ll) reported higher X-ray luminosi- 
ties f or the non-accre t ing W TTS stars than for CTTS. In the 
ONC. iFeigelson et al. (2002) concluded from Chandra obser- 
vations that the presence of circumst ellar disks has no influ - 
ence on the X-ray emission, whereas [Flaccomio et al. ( 2003a ). 
in another Chandra study of the ONC, found Lx and Lx/L* to 
be enhanced in WTTS when compared to CTT S. From the re- 



cent C handra Orion Ultradeep Project (COUP). lPreibisch et al 



(2005) reported the X-ray emission of WTTS to be consistent 
with the X-ray emission of active Main Sequence (MS) stars, 
while it is suppressed in CTTS. However, in all these studies 
the X-ray emission mechanism is consistent with a scaled-up 
version of a solar corona. 

X-ray emission during accretion outbursts has been ob- 
served in VI 647 Or i dKastner et all 12004 iGrosso et all 120051: 
Kastner et all 120061) and in VI 1 18 Ori dAudard et all l2005b . 
The X-ray luminosity increased by a factor of 50 during the 
outburst in VI 647 Ori, and the spectrum hardened. On the other 
hand, the X-ray luminosity of VI 1 18 Ori remained at the same 



level as during the pre-outburst phase, while the spectrum be- 
came softer. 

Possible signs of accretion-induced X-ray emission are 
revealed in a few high-resolution spectra of CTTS. High 
electron densities were measu red in the spectra of TW 



Hva (IKastner et all 120021 Istelzer & Schmiti 12004 . BP Tau 
dSchmitt et all \20oA iRobrade & Schmittl l2006h ."and V4046 
Sgr dGiinfher et all 120061) . and were interpreted as indications 
of X-ray production in accretion shocks. Other spectroscopic 
features that also suggested an accretion shock scenario are the 
low electron temperature dominating the plasma in TW Hya (a 
few MK, as exp ected from shock-induced heating) and abun- 
dance anomalies. Istelzer & Schmittl d2004 interpreted the high 
Ne/Fe abundance ratio as being due to Fe depletion by cond en- 
sation into grain in the accretion disk. iDrake et al. (2005) re- 
ported a substantially larger Ne/O ratio in the spectrum of TW 
Hya than in the spectra of the other studied stars, and they pro- 
posed to use this ratio as a diagnostic for metal depletion in the 
circumstellar disk of accreting stars. 

Work on high-reso lution X-ray s pectros copy was subse- 
quently extended by iTelleschi et all d2006al) to a sample of 
9 pre-main sequence stars with different accretion properties. 
The main result of that work is the identification of an excess 
of cool plasma measured in the accreting stars, when compared 
to WTTS. The origin of this soft excess is unclear. Further ev- 
idence for a strong soft excess in the CT TS is revealed in th e 
extraordinary X-ray spectrum of T Tau dGiidel et all |2006cJ). 
In this case, however, the electron density (derived from spec- 
tral lines formed at low temperatures) is low, n e < 10 10 cirT 3 . 
The density, in case of accretion shocks, can be estimated us- 
ing the strong shock condition «2 = 4«i, where n\ and «2 
are the pre-shock and post-shock densities, respectively. The 
density n\ can be derived from the accretion mass rate and 
the accreting area on the stellar surface: M w 4nR 2 fvffn e m p , 
where / is the surface filling factor of the accretion flow, and 
Vff = (2GM/R)^ 2 is the free-fall velocity. Usin g the accretion 
rate M « (3 - 6) x IP' 8 M G yr 1 for T Tau dWhite&Ghez . 
l200lHCalvet etalll2004 we obtain n 2 = (1.1 - 2.2) x 10 11 // 
dGiidel et all 2006c ). Even in the extreme case that / = 10%, 
we expect a density > 10 12 cm -3 , i.e. orders of magnitude 
higher than the measured value. 

The aim of the present paper is to study the role of accre- 
tion in the overall X-ray properties of pre-main sequence stars 
in the Taurus-Auriga Molecular Cloud, by coherently compar- 
ing samples of CTTS and WTTS. Our analysis is complemen- 
tary to the COUP survey work, and we will present our re- 
sults along largely similar lines (see IPreibisch et alJl2005l for 
COUP). Indeed one of the main purposes of the present work 
is a qualitative comparison of the Taurus results with those ob- 
tained from the Orion sample. We do not, however, present is- 
sues related to rotation; rotation-activi ty relations will be sepa- 
rately discussed in a dedicated paper (Brigg s etalll2006h . 

The paper is structured as follows. In Sect. [2] we describe 
the stellar sample used in this work, and in[3]we summarize the 
relevant steps of the data reduction. We present our results in 
Sect. HI and discuss them in Sect. [5] We summarize our results 
and conclude in Sect. [6] 
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2. Studying X-rays in the Taurus Molecular Cloud 

2.1. The Taurus Molecular Cloud 

We will address questions on X-ray production in accret- 
ing and non-accreting T Tauri stars using data from the 
XMM- Newton Extended Su rvey of the Taurus Molecular Cloud 
(XEST, iGiidel et al.ll2006al) . 

The Taurus Molecular Cloud (TMC) varies in significant 
ways from the Orion Nebula Cluster and makes our study an 
important complement to the COUP survey. The TMC has, 
as the nearest large , star-forming region (distance » 140 pc, 
lLoinard et alj|2005l) . played a fundamental role in our under- 
standing of low-mass star formation. It features several loosely 
associated but otherwise rather isolated molecular cores, each 
of which produces one or only a few low-mass stars, different 
from the much denser cores in p Oph or in Orion. TMC shows 



a low stellar density of only 1-10 stars pc 2 (e.g JGomez et al 



19931) . In contrast to the very dense environment in the Orion 



Nebula Cluster, strong mutual influence due to outflows, jets, or 
gravitational effects is therefore minimized. Also, strong winds 
and UV radiation fields of OB stars are present in Orion but 
absent in the TMC. 

The TMC has provided the best-characterized sample of 
CTTS and WTTS, many of which have b een subject to detailed 
studie s; see, e.g., the seminal work by iKenvon & Hartmann 
d 1995b that concerns, among other things, the evolutionary 



history of T Tau stars and their disk+envelope environment, 
mostly based on optical and infrared observations. 

It is therefore little surprising that comprehensive 
X-ray studies of selected objects as well as surveys 
have been performed with several previous X-ray satel- 



Feiselson et al. dl987l). Walter etal 


dl988 


). Bouvier!dl990l). 


Strom et al.ldl990h.lStrom & Stromll 


1994. 


Damiani & Micela 



(1995), Damiani etal. (1995J, iNeuhauser et all (119951) . and 
Stelzer & Neuhauser (2001). Issues we are studying in our 



paper have variously been stud ied i n these surveys before, 
although, as argued by IGiidel et alJ (I2006al) and below, the 
present XEST project is more sensitive and provides us with 
a near-complete sample of X-ray detected TTS in the surveyed 
area, thus minimizing selection and detection bias. 

2.2. The XEST sample of T Tau stars 

The XEST project is an X-ray study of the most populated re- 
gions (comprising an area of « 5 square degrees) of the Taurus 
Molecular Cloud. The survey consists of 28 XMM-Newton ex- 
posures. The 19 initial observations of the project (of approxi- 
mately 30 ks duration each, see Table 1 of IGiidel et aD Eo06a) 



were complemented by 9 exposures from other projects or from 
the archive. Also, 6 Chandra observations have been used in 
XEST, to add information on a few sour ces not detected wi th 
XMM-Newton, or binary information (see lGiidel et alj |2006a). 

To distinguish between CTTS and WT TS we use the clas - 
sification given in col. 10 of Table 1 1 of IGiidel et alJ (I2006al) . 
This classification is substantially based on the equivalent 
width of the Ha line (EW[Ha-]). For spectral types G and K, 



stars with EW(Hff) > 5 A are defined as CTTS, while other 
stars are defined as WTTS. For early-M spectral types, the 
boundary between CTTS and WTTS was set at EW(Ha) = 
10 A and for mid-M spectral types at EW(Ha) = 20 A. Stars 
with late-M spectral type are mainly Brown Dwarfs (BDs), and 
given their low optical continuum, a clear accretion criterion 
is difficult to provide. For this reason, BDs were treated as a 
class of their own and are not used in our comparison studies 
of accretors vs. non-accretors (but were included i n the "tota l" 



samples wh en appropriate ) . For further details, see Giideletal 



d2006al) and lGrosso et all d2006l) . YSO IR types were used to 
classify borderline cases and protostars. In summary, protostars 
have been classified as type or 1 (Class and I, respectively), 
CTTS are type 2 objects, WTTS are type 3 objects, and BDs are 
classified as type 4. Type 5 is assigned to Herbig Ae/Be stars, 
while stars with uncertain classification are assigned to type 9. 
We will use these designations in our illustrations below. 

We emphasize the near-com pleteness of XEST with re- 
gard to X-ray detections of TTS. Gtid el et al.l d2006al) provide 
the detection statistics (their Table 12): A total of 126 out of 
the 159 TMC members surveyed with XMM-Newton were de- 
tected in X-rays. Among these are 55 detected CTTS and 49 
detected WTTS (out of the 65 and 50 surveyed targets), corre- 
sponding to a detection fraction of 85% and 98%, respectively. 
Almost all objects have been found comfortably above the ap- 
proximate detection limit of Lx ~ 10 28 erg s , indicating that 
TTS generally emit at levels between 10 29 - 10 31 erg s _1 , ex- 
ceptions being lowest-mass stars and brown dwarfs. Most of 
the non-detected objects have been recognized as stars that are 
strongly a bsorbed (e.g. , by thei r own disks) or as stars of very 



low mass dGiidel et al. . 2006a ). XEST is the first X-ray sur- 



vey of TMC that reaches completeness fractions near unity, 
and therefore minimizes detection bias and unknown effects of 
upper limits to correlation studies as performed here. It pro- 
vides, in_fhisjegard, an ide al comparison with the COUP re- 



2005). 



suits ( Prei bisch et al 

A few sources were excluded from consideration in the 
present work. These are the four stars that show composite X- 
ray spectra possibly originating from two differe nt sources (DG 
Tau A, GV Tau, DP Tau, and CW Tau; IGiidel et al.l[2006bh . 
and three stars which show a decreasing light curve throughout 
the observation (DH Tau, FS Tau AC, and V830 Tau); these 
light curves probably describe the late phases of large flares. 
Further, the deeply embedded protostar L1551 IRS5, which 
shows ligh tly absorbed X-ray emission that may be attributed 
to the jets dFavata et all 120021 iBally et all 120031) . was also ex- 
cluded, and so were the two Herbig stars (AB Aur, V892 Tau). 
In some correlation studies, we do conside r objects for which 
upper limits to Lx have been estimated in IGiidel et al.l d2006al) . 
but will not consider non-detections without such estimate (as, 
for example, if the absorption is unknown). 

Our final, basic sample of TTS then consists of 56 CTTS 
and 49 WTTS. Among the X-ray detections, there are also 8 
protostars, 8 BDs, 2 Herbig stars and 4 stars with uncertain 
classification. Smaller subsamples may be used if parameters 
of interest were not available. 

When L, is involved in a correlation, we excluded 
all stars that are apparently located below the Zero- Age 
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Main Sequ ence (ZAMS) in th e Hertzsprung-Russel diagram 
(Fig. 10 of lGudeletaf]l2006ah . These stars are HH 30, IRAS 
S04301+261, Haro 6-5 B, HBC 353, and HBC 352. Their lo- 
cation in the HRD is likely to be due to inaccurate photometry. 

Many of the stellar counterparts to our X-ray sources 
are unresolved binaries or multiples. In total, 45 out of the 
159 stellar systems surveyed by XMM-Newton are multiple 
dGudel et alll2006al) . f - as we will find below in general, and 
as has been reported in earlier studies of T Tauri stars (e.g., 
Prei bisch et all 120051) - Lx scales with L», then this also holds 



for the sum of the Lx with respect to the sum of the component 
L„. Binarity does therefore not influence comparisons between 
Lx and L*. When correlating Lx with stellar mass, we will find 
that more massive stars are in general brighter. In the case of 
binaries, the more massive component (usually the more lumi- 
nous "primary" star) will thus dominate the X-ray emission. 
We have used the primary mass for the stellar systems if avail- 
able; we therefore expect the influence of the companions on 
our correlations to be small. We will also present tests with the 
subsample of single stars below. 

3. Data reduction and analysis 

The XEST survey is principally based on CCD camera expo- 
sures, but is complem ented with high resolut ion grating spectra 
for a few bright stars dTelleschi et allr2006ah . and with Optical 



Monitor observations ( Aud ard et al.L 12006). The three EPICs 



onboard XMM-Newton are CCD-based X-ray cameras that col- 
lect photons from the three telescope s. Two EPIC detectors are 
of the MPS type dTurner et allEoOlh and one is of the PN type 
(Strii der et all l2Q0lh . They are sensitive in the energy range of 

0. 15-15 keV with a spectral resolving power of E/AE = 20-50. 

The data were reduced using the Science Analysis System 
(SAS) version 6.1. A detailed des cription of all data reduction 
procedures is given in Sect. 4 of Gudel et al. ( 2006ah . 

Source and background spectra have been obtained for each 
instrument using data during the Good Time Intervals (GTIs, 

1. e., intervals that do not include flaring background). Further, 
time intervals with obvious strong stellar flares were also ex- 
cluded from the spectra in order to avoid bias of our results by 
episodically heated very hot plasma. 

One PN spectrum usually provides more counts than the 
two MOS spectra together. We therefore only used the PN data 
for the spectral analysis, except for the sources for which PN 
data were not available (e.g., because the PN was not opera- 
tional, or the sources fell into a PN CCD gap). 

The spectral fits were performed using two different ap- 
proaches in the full energy band. First, we have fitted the spec- 
tra using a conventional one- or two-component spectral model 
(1-T and 2-T), both components being subject to a common 
photoelectric absorption. In this approach, the hydrogen col- 
umn density, A^h, two temperatures (Ti t} and two em ission 
measures (EMi^) are fitted in XSPEC dArnaudl Il996h using 
the vapec thermal collisional-ionization equilibrium model. 

In the second approach, the spectra were fitted with 
a model consisting of a continuous emission measure dis- 
tribution (EMD) as found for pre-main sequence and ac- 
tive ZAMS stars dTelleschi et all 120051: lArgiroffi et all [2004; 



Garcia- Alvarez et all |2005e IScelsi et all 120051) . The model 
consists of a grid of 20 thermal components binned to intervals 
of J log T = 0.1 from logT = 6 to log T =7.9, arranged such 
that they form an EMD with a peak at a temperature To and two 
power-laws toward lower and higher temperatures with power- 
law indices a and j3, respectively. Given the poor sensitivity of 
CCD spectra at low temperatures, a was kept fixed at 2, con 



sistent with values found in previous studies dTelleschi et al 
20051 lArgiroffi et all |2004 . while we let /? free to vary (be- 
tween -3 < jS < 1). The absorbing hydrogen column density 
A^h was also fitted to the data. The abundances were fixed at 
values typical for pre-ma in sequence stars or very active zero- 
age main-sequence stars ((Telleschi et all 120051: [ Argiroffi et all 



2004; Garcia- A 



ther details, see 



varez et alll2005blScelsi et all l2005H'l Forfur- 
Giidel et all d2006al)~ 



For each star and each model we computed the average 
temperature (r av ) as the logarithmic average of all temperatures 
used in the fit, applying the emission measures as weights. The 
X-ray luminosity (Lx) was computed in the energy range 0.3- 
10 keV from the best-fit model assuming a distance of 140 pc. 

Of the 126 members detected in XEST, 22 were detected 
in two different exposures. In those cases, two separate spectral 
fits were made. For correlations of X-ray parameters with stel- 
lar properties, we used logarithmic averages of the results from 
the two fits. On the other hand, if we correlate X-ray properties 
with each other, we treat the two spectral fit results from the 
same source as different entries. 

Results from the spectral fits are given in Tabl e 5 (for the 



EMD fits) and Table 6 (for the 1-T and 2-T fits) of Gu del et al 



(2006a). We use the results from the EMD interpretation to per- 
form statistical correlations below. 



4. Results 

Motivated by results from previo us X-ray studies and in partic- 
ular guided by the COUP work dPreibisch et all l2005). we now 
seek systematics in the X-ray emission by correlating X-ray 
parameters first with fundamental stellar parameters, and then 
also seeking correlations among the X-ray parameters them- 
selves. We will consider the fundamental stellar properties of 
mass and bolometric luminosity, accretion rate, and age, but 



we w ill not discuss rotation properties here (see iBriggs et al 



120061 for a detailed study). The basic X-ray properties used for 
our correlations are the X-ray luminosity Lx in the 0.3-10 keV 
band and the average electron temperature T av . One of the main 
goals of this section is to seek differences between CTTS and 
WTTS. We will compare our findings with those of COUP and 
some other previous work in Sect. [5] 



1 The abundance values used are, with respect to the solar 
photospheric abundances of Ande rs & Grevess J Jl989h : C=0.45, 
N=0.788, 0=0.426, Ne=0.832, Mg=0.263, Al=0.5, Si=0.309, 
S=0.417, Ar=0.55, Ca=0.195, Fe=0.195, Ni=0.195 
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4.1. Correlations between X-rays and stellar 
parameters 

4.1 .1 . Correlation with mass 

In Fig.Q]we plot the X-ray luminosity, Lx (in erg s _1 ), as a func- 
tion of the s tellar mass (M, in units of the solar mass, M , from 
Table 10 in lGiidel et alj|2006al) . In the left panel we show the 
relation for all types of objects in our sample. Different sym- 
bols are used to mark different object types (see panel in fig- 
ure). Upper limits for non-detections are marked with arrows. 
In the middle and right panels the same relation is shown sep- 
arately for CTTS and WTTS, respectively. A clear correlation 
is found between the two parameters in all three plots, in the 
sense that Lx increases with mass. The correlation coefficients 
are 0.79 for the whole sample (99 entries), 0.74 for the CTTS 
(45 entries), and 0.84 for the WTTS (43 entries). We com- 
puted the sig nificance of the corr elation using correlation tests 
in ASURV dLa Valley etaTlll992t specifically, the Cox hazard 
model, Kendall's tau, and Spearman's rho have been used) and 
found a probability of < 0.01% that the the parameters are un- 
correlated in each of the three cases. As for all subsequent sta- 
tistical correlation studies, we summarize these parameters in 
Table □ 

We computed linear regression functions for the loga- 
rithms of the two parameters, of the form logy - a + b\ogx, 
using the parametric estimation maximization (EM) algo- 
rith m in ASURV, which implements the methods presented by 
Isobe et alJ (119861) . We find the regression functions logLx = 
(1.69 + 0.1 1) logM + (30.33 + 0.06) for the full stellar sample, 
logLx = (1-70 ± 0.20) log M + (30.13 + 0.09) for the CTTS, 
and logLx = (1-78 ± 0.17) log M + (30.57 + 0.09) for the 
WTTS. The regression parameters are also listed in Table Q] 
as f or all subsequent regr ession analyses. In the ONC sam- 
ple, IPreibisch et al.1 (120051) found the similar linear regression 
logLx = (1.44 + 0.10) logM + (30.37 + 0.06) for all stars with 
masses < 2M using the same algorithm. 

The EM algorithm is an ordinary least-square (OLS) re- 
gression of the dependent variable y (Lx in this case) against 



the independent variable x (M). When using this method, we 
assume that Lx is functionally dependent on the given mass 
dlsobe et al. , 1990). However, the M values are also uncertain, 
and assuming a functional dependence a priori may not be cor- 
rect. We therefore also comp uted the lin e ar reg ression using 
the bisector OLS method after llsobe et al.l dl990 ). which treats 
the variables symmetrically. In this case, we find logL x = 
(1.91 ± 0.11) log M + (30.44 + 0.05) for all stars together, 
logL x = (1.98 + 0.20)logM + (30.24 + 0.06) for the CTTS, and 
logLx = (2.08+0.17) log M+ (30.69 +0.07) for the WTTS (see 
Table [U. The slopes for the bisector OLS are slightly steeper 
than in the EM algorithm. However, the values for CTTS and 
WTTS agree within one sigma, and we caution that the upper 
limits for non-detections are not taken into account in the bi- 
sector linear regression method. 

We verified this trend for the subsample of stars that have 
not been recognized as multiples. We find the regression lines 
logLx = (1-72 + 0.12)logM + (30.39 + 0.07) using the EM 
algorithm, and logL x = (1.85 ± 0.11)logM + (30.48 ± 0.07) 
using the bisector algorithm. These results are fully consistent 
with the results for the total sample. 

Differences are present between the CTTS and WTTS stel- 
lar samples. While the slopes found in the correlations for 
CTTS and WTTS are consistent within lcr, the intercept of 
WTTS at lM (log M/M Q = 0) is * 0.45 dex larger than the 
intercept for CTTS. This lets us anticipate a larger average 
Lx in WTTS. The correlation is better determined for WTTS, 
as judged from a slightly higher correlation coefficient and a 
smaller error in the slope. Furthermore, the standard deviation, 
cr, of the points with respect to the regression function from the 
EM algorithm is slightly larger for CTTS (0.45) than for WTTS 
(0.38). 

4.1 .2. Evolution of X-ray emission 

Here, we discuss the evolution of the X-ray emission with age. 
Among main-sequence (MS) stars, Lx is correlated with rota- 
tion and anti-correlated with age. The common explanation is 
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Table 1. Summary of results found for the different correlations. In the third column, n is the number of stars used in the statistic. 
P is the probability that the parameters are uncorrelated (computed with ASURV), C is the correlation coefficient and cr is the 
standard deviation from the EM algorithm. The intercept of the linear regression is a, and b is the slope. Errors are 1-sigma rms 
values for the respective variables. 
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Fig. 2. Evolution of the normalized X-ray emission. Lx has 
been normalized with the predicted values from the Lx-mass 
relation (see text for details). The linear regression computed 
with the EM algorithm is plotted (solid line) together with er- 
rors in the slopes (dashed line). Symbols mark different types 
of stars. 



that magnetic activity is directly related to the stellar rotation, 
and the latter decays with age because of magnetic braking. 
However, TMC PMS stars do not sh ow the relation betw een 
Lx and rotation observed for MS stars (Bri ggs etalll2006h .On 
the other hand, Lx decreases during the evolution of pre-main 
sequence stars, provided that a common X-ray saturation law 
applies (see below), because L* decreases along the Hayashi 
track. 

We have found (see Sect. |4. 1 .Tb that L x shows a strong cor- 
relation with mass. In order to avoid an interrelationship be- 
tween the correlations, we normalize the measured Lx with 



Lx(M) predicted by the correlation with mass (Lx(M) = 
jq30.33^i .69 er g s -i) an( j mu itiply with the Lx expected for 

a 1 M Q star (10 3033 erg s _1 ). We designate this quantity by 
L X (M = 1M ). In Fig. |2] we plot L X (M = 1M ) as a function 
of age. A slight decline in L x is found between 0. 1 and 10 Myr. 
The correlation coefficient is C = -0.31 for 93 entries. The 
tests in ASURV give probabilities between 0.1% and 0.38% 
that age and Lx(M = 1M ) are uncorrelated. We have com- 
puted a linear regression with the EM algorithm in ASURV and 
find log(L x /L x (M)) = (-0.36 + 0.11) log (age) + (30.45 +0.06) 
(where age is in Myr). Further, we have tested the linear regres- 
sion and the correlation probability when we neglect the two 
youngest stars (V410 X4 and LkHa 358) and found the linear 
regression to be consistent within error bars with the above re- 
lation, with a probability of P < 1 % for no correlation. Further, 
as a test, we have computed the linear regression using the bi- 
sector algorithm. We find a much steeper slope of -1.02 + 0.07 
(Table [TJ, indicating that the linear regression is nevertheless 
only marginal, and the scatter is dominated by other contribu- 
tion^ 

Preibisch & Feigelsonl ( 2005) reported correlations consis- 



tent with ours, applying the EM algorithm to the ONC data. 
They correlated Lx with age in mass-stratified subsamples of 
the surveyed stars and found Lx to decrease with age with 
slopes ranging from -0.2 to -0.5, i.e. fully consistent with the 
slope found in Fig. [2] (mass-stratified analysis for our sample 
also indicates decreasing Lx in some mass bins but not in oth- 
ers; our statistics are too poor for this purpose). 

4.1 .3. Mass accretion rates 

We now directly compare the X-ray parameters derived from 
our spectral fits (Lx, Lx/L t , and L av ) with the previously deter- 
mined mass accret ion rates (M, in M y r' 1 ). We use M listed in 
the XEST catalog (Gudeletal., 2006a, and references therein). 
Accretion rates may be variable, and various methods for their 
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Fig. 3. Residual X-ray luminosity for CTTS (after normaliza- 
tion with the M - Lx and M - M relation) as a function of the 
mass accretion rate. Regression lines obtained using the EM al- 
gorithm (red) and the bisector algorithm (blue) are plotted with 
their respective errors in the slope (dashed lines). The regres- 
sion lines are computed only using the stars plotted with black 
circles, while the two stars with M smaller than 10~ 10 M o yr~' 
were ignored (see text for more details). 

determination may produce somewhat different results. If dif- 
ferent values were found for a given star in the literature, the 
range of M is marked by a horizontal line in our figures. 

When comparing Lx with M, some caution is in order. In 
Sect. 14. 1.1 1 we have shown that a tight relation exists between 
L x and the stellar mass. Further, a clear relation between M and 
the stellar mass has been found for class I I objects in the liter - 



ature (e.g., iMuzeroUe et all [2003, 2005; ICalvet et all 120041) . 



Combining the Lx-mass and M-mass relations, we expect Lx 
to correlate with M as well. However, here we are interested 
in testing if an intrinsic relation between the latter two parame- 
ters exists that is no t a consequence of the two form er relations 



Calvet et al. (2004) have used evolutionary tracks of lSiess et al 
(|200pJ) (consistent with the XEST survey) to find a relation 
MocM 195 in the mass rang e between 0.02 and 3 M B . Similarly, 
Muzerolle etail d2003l) and lMuzerolle et all d2005l) found M oc 
M 2 and M oc M 2A respectively, using different evolutionary 
tracks. We therefore adopt the relation logM « 2 logM - 7.5. 
Further, we use logL x = 1.69 logM + 30.33 (Sect. |47T7TT >. and 
we then compute the expected Lx for each M value, namely 
log L X (M) = 0.85 logM + 36.67. 

In Fig. [3] we plot the ratio of Lx/Lx(M) as a function of 
M for class II objects. We will expect that the values scatter 
around a constant if this ratio were determined only by the 
M - Lx and M - M relations. We find a very large scatter for 
any given M (2-3 orders of magnitude) but, using a regres- 
sion analysis, a tendency for weak accretors to show higher 
Lx, compared to strong accretors. However, if we exclude the 
two stars with the smallest accretion rates (plotted with blue 
symbols in Fig. |3), the correlation is less clear. In this case, 
the correlation coefficient is C — -0.47 for 37 data points. 
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Fig. 4. Upper panel: Fractional X-ray luminosity vs. M. Lower 
panel: T- dv vs. M. Different symbol sizes represent different ob- 
ject types as defined in the figures. Arrows represent upper lim- 
its for the accretion rates. 



Nevertheless, the probability, computed in ASURV for the EM 
algorithm, that no correlation is present is only P = 0.1% - 

0. 5%. We have computed the linear regression using differ- 
ent methods. With the EM algorithm we find log Lx/Lx(M) = 
(-0.48 ±0.15) logM -(4.05 + 1.19) (Tabled). Using the bisec- 
tor algorithm, however, the slope is found to be -1.02 + 0.14, 

1. e., more than 3 sigma steeper than the slope found with the 
EM algorithm. The entries for two stars with low M (plotted in 
blue) are consistent with the linear regression found with the 
EM algorithm. We conclude that the two parameters are not 
evenly distributed, but that a linear regression of the logarith- 
mic values cannot clearly be claimed. 

In Fig.[4]we plot Lx/L t and T av as a function of the accre- 
tion rate. Arrows represent upper limits for M. In both cases no 
correlation is evident. 
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L.Ao 

Fig. 5. Lx as a function of the L». Left: for all stars; the different symbols describe different classes of stars, while arrows are 
upper limits for non-detections. Middle: same for CTTS (type 2). Right: same for WTTS (type 3). The horizontal bars show the 
ranges of literature values for L» . 



4.1.4. Correlation with bolometric luminosity 

In Fig. [5] we plot L x as a function of t he stellar bolometric lu- 
minosity L, (from iGiidel et al. 2006a and references therein). 



The lines corresponding to L x /L* = 10 , 10 , and 10 are 
also shown. In the left panel, all stars are plotted, with differ- 
ent symbols for each stellar class as described in the figure. We 
again excluded from the plot the stars mentioned in Sect. [2] 
Upper limits for non-detections are marked with arrows. By 
far most of the stars are located between L x /L» = 10~ 3 and 
Lx/L* = 10 -4 . In the middle and right panels, we present CTTS 
and WTTS separately. 

The correlation coefficients are 0.83 for the full stellar sam- 
ple (108 entries), and 0.84 (48 entries) and 0.85 (44 entries) 
for CTTS and WTTS, respectively. Probabilities for the ab- 
sence of a correlation are very small, P < 0.01%. We computed 
linear regression lines with the EM algorithm in ASURV. For 
the full sample, we find logL x = (1.05 ± 0.06) log L*/L + 
(30.00 + 0.05), whereas for CTTS and WTTS, logL x = 
(1.16 ±0.09) log L»/L Q + (29.83 + 0.06) andlogL x = (1.06 ± 
0. 1 0) log L t /Lq + (30.22 + 0.08), respectively. Fig. EJ; shows one 
WTTS at LJL e * 0.01 with a rather high L x ~ 5 x 10 29 erg s" 1 
(KPNO-Tau 8 = XEST-09-022). Not considering this object, 
the slope of the regression slightly steepens to 1.17 + 0.09, 
which is only marginally different from the slope based on all 
WTTS. The standard deviation at the same time marginally de- 
creases from 0.41 to 0.36. 

Again, we also computed a linear regression using the 
bisector OLS algorithm that treats both L, and L x as inde- 
pendent variables. The slopes are very similar to the those 
reported above (Table [T). The important distinction between 
CTTS and WTTS is that the latter clearly tend to be located at 
higher L x /L» (see also below): at L*/L = 1 the CTTS show 
an average logL x = 29.83 [erg s _1 ], while for the WTTS, 
logL x = 30.22 [ergs- 1 ]. 

The regressions are thus compatible with a linear relation 
between L x and L«, and therefore L x /L» is, on average for a 



given L„, a constant between 10~ 4 and 10~ 3 regardless of L t . 
This is reminiscent of the situation among very active, rapidly 
rotating main-sequence or evolved subgiant stars that saturate 
at fractional X-ray luminosities of the same order, provided 
they rotate sufficiently rapidly. We thus find that the majority of 
our TTS are in a saturated state. A consequence of this would 
be that ro tation no longer cont rols the X-ray output, as sug- 
gested by Preibisc h et al. (2005) f or the Orion s ample . This is 
discussed for the XEST sample by Briggs et al.1 (2006). Below, 
we will specifically study whether the Lx/L* relation is differ- 
ent for CTTS and for WTTS. 

The correlation found for the full stellar sample is com- 
patible with the relation found in Orion: log L x = (1 .04 + 



0.06)log(L»/Lo) + (30.00 + 0.04) dPreibisch et al.. 2005). The 



slope found for the WTTS in the ONC is also consistent with 
our results within the error bars. For CTTS, on the other hand, 
Preibisch et al. d2005l) bund a very large scatter in the correla- 
tion. This is not observed in our XEST sample; we rather see 
similar scatter for CTTS and WTTS, as demonstrated by the 
similar correlation coefficients, the similar errors in the s lope, 
and the similar standard deviations. IPreibisch et al.l (12005) sug- 
gested that strong accretion could lead to larger errors in the 
determination of stellar luminosity and the effective tempera- 
ture. 

4.1 .5. Fractional X-ray luminosity L x /X* 

In Fig. [6] we plot the histogram for the distribution of 
log(Lx/X») for CTTS (grey) and WTTS (white). The two 
populations are different, with WTTS having a larger mean 
log(L x /L„). We fitted each of the two histograms with a 
Gaussian function and computed the mean and its errors. For 
CTTS, we find <log(L x /L„)> = -3.73 + 0.05, while for WTTS, 
<log(L x /L*)> = -3.39 + 0.06. 

A more rigorous test is based on the Kaplan-Meier estima- 
tors as comp uted in ASURV, wh i ch im plements the methods 
presented bv lFeigelson & Nelson! dl985l) . This method also ac- 
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Fig. 8. Cumulative distributions of log(L x /L,,) for WTTS (dotted) and CTTS (solid) for masses smaller than 0.3 M Q (left), masses 
between 0.3 and 0.7 M Q (middle), and masses larger than 0.7 M Q (right). 
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Fig. 6. Distributions of log(Lx/L») for CTTS (grey histogram) 
and WTTS (white histogram). Gaussians that fit the distribu- 
tions are also plotted and their peaks are marked with vertical 
lines. 



counts for the upper limits in L x for the non-detections. The 
results are plotted in Fig. [7] The solid line represents the CTTS, 
the dotted line the WTTS. The WTTS distribution is shifted to- 
ward larger log(L x /L„,) compared to the CTTS distribution by a 
factor of approximately 2. We find (log(L x /L*)) = -3.72+0.06 
and <log(L x /L,)> = -3.36 + 0.07 for CTTS and WTTS, re- 
spectively, in full agreement with the Gaussian fit. Judged from 
a two-sample test based on the Wilcoxon test and logrank test 
in ASURV, the probability that the two distributions are ob- 
tained from the same parent population is very low, namely 
P = 0.01%-0.03%. 

Again, we test this result using the subsample of stars that 
have not been recognized as multiples. The subsample consists 
on 29 CTTS (4 of which have upper limits) and 33 WTTS (with 
no upper limits). We find a probability of P = 0.05%-0.07% 
that the distributions arise from the same parent population, 



Fig. 7. Cumulative distribution of log(L x /L.) for CTTS (solid) 
and for WTTS (dotted). 



and (log(L x /L«)>cTTS = -3.83 + 0.06 and (log(L x /L«)>WTTS = 
-3.40 + 0.08. These results are consistent with the results found 
in the full sample. We can therefore conclude that multiplicity 
does not influence our results. 

In Fig. [8] we plot the Kaplan-Meier estimator for the dis- 
tribution of log{Lx/L„) in three different mass ranges: for stars 
with masses smaller than 0.3 M Q , between 0.3 and 0.7 M , and 
larger than 0.7 M B . For the latter two mass bins, the distribu- 
tions belong to two different parent populations at the > 94% 
level. For lower masses (M < 0.3M G ), we find the CTTS 
and WTTS distributions to be similar, but we still find larger 
log(L x /L*) for WTTS than than for CTTS. The difference be- 
tween the two populations is not significant at the « 10% level 
possibly because of the small size of the stellar sample in this 
mass range. 

Our results can be compared with t he distributions found i n 
the COUP survey, shown in Fig. 16 of Preibisc h et al. (2005). 
In the latter figure, the stars are classified according to the 8542 
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Fig. 10. X-ray Luminosity Function (XLF) for CTTS (solid) Fig. 11. Stellar bolometric luminosity function for CTTS 
and WTTS (dotted). (solid) and WTTS (dotted). 



A Ca II line, which is an indicator of disk accretion, similar to 
the EW(Ho') used in our work. In Orion, a substantial difference 
has been found between the distributions of accreting and non- 
accreting stars in the mass ranges 0.2-0.3 M and 0.3-0.5 M . 
However, for 0.5-1 M , the two distributions appeared to be 
compatible, in contrast to our findings that show fainter CTTS 
consistently in all mass ranges. 

4.2. Correlations between X-ray parameters 

4.2.1 . The X-ray luminosity function 

In Fig.[l0]we display the X-ray luminosity function (XLF) for 
WTTS and CTTS for our Taurus sample. The XLF has again 
been calculated using the Kaplan-Meier estimator in ASURV, 
so that the few upper limits have also been considered. The 
total number of sources used was 105, 56 of them being CTTS 
(including 6 upper limits) and 49 WTTS (including 1 upper 
limit). The WTTS are again more luminous than CTTS by a 
factor of about 2 (with mean values (logLx)c = 29.51 and 
(logLx)w = 29.80). The probability that the two distributions 
arise from the same parent population is 7%-10%, computed 
using ASURV as described in Sect. 14.1.51 

If we restrict the stellar sample to stars with no recog- 
nized multiplicity, we obtain average X-ray luminosities of 
(logLx)c = 29.38 and (logLx)w = 29.65, for samples con- 
sisting of 32 CTTS (5 upper limits) and 36 WTTS (1 upper 
limit). The difference between the two stellar samples is 0.3 
dex (i.e., a factor of two), similar to what we found for the full 
sample. However, the two-sample tests give a larger probabil- 
ity (P = 12%— 33%) that the two stellar groups arise from the 
same parent population. Among the multiple sources, 24 are 
CTTS (with 1 upper limit), but only 13 are WTTS (with no up- 
per limit). By adding the multiples to the sample of single stars, 
we expect that the distributions slightly shift toward larger Lx, 
and because there are significantly more multiple CTTS, the 
CTTS distribution of the total sample should be more similar 



to the WTTS total distribution, but the opposite trend is seen. 
We conclude that the trends are grossly the same for the total 
sample and the single-star subsamples, the larger probability 
being due to significantly smaller samples that are compared. 
Overall, thus, CTTS are recognized as being X-ray deficient 
when compared to WTTS. 

Fig. [9] shows the X-ray luminosity function for the same 
three mass ranges as used in Sect. 14.1.51 (M < O.3M in the 
left panel, O.3M < M < O.7M in the middle panel, and 
M > O.7 in the right panel). Again, we find the largest differ- 
ence between CTTS and WTTS for the two higher-mass bins, 
with probabilities of only 2%-4% that the distributions belong 
to the same parent population. The probability is substantially 
larger for M < O.3M (29%-32%), but the statistics are also 
considerably poorer. 

Considering the difference in the XLFs of CTTS and WTTS 
alone, a possible cause could be that the bolometric luminos- 
ity function of CTTS would indicate lower luminosities L„ 
than for WTTS, which would result in lower average Lx pro- 
vided that Lx/L„ « constant, i.e., that saturation applies for all 
stars. In Fig.QT| we plot the distributions of L* for WTTS and 
CTTS. In fact, the CTTS are found to be slightly more lumi- 
nous than the WTTS. We find <logL*> c = 33.35 + 0.08 and 
(logL,) w = 33.19 + 0.09. The probability that the distribu- 
tions arise from the same parent population is 15%— 21 %, i.e. 
making the difference marginal. We conclude that because Lx 
is linearly correlated with L, (Fig. |5), the difference in Lx for 
the two samples is intrinsic, which is of course a reconfirma- 
tion of our previous finding that the distributions of Lx/L, also 
indicate lower activity for CTTS compared to WTTS. 

4.2.2. Absorption 

In Fig. Q~2] we plot the distribution of Nh for accreting (solid 
line) and non-accreting stars (dotted line) calculated using the 
Kaplan-Meier estimator in ASURV. The logarithmic average 
of N H for CTTS is (Nn) c = 4.2 x 10 21 cirT 2 and is more than 
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Fig. 9. X-ray Luminosity Function (XLF) for CTTS (solid) and WTTS (dotted) for for masses smaller than 0.3 M (left), masses 
between 0.3 and 0.7 M (middle), and masses larger than 0.7 M (right). 




Fig. 12. Cumulative distribution of Nh for CTTS (solid line) 
and WTTS (dotted line). 

a factor of two larger than the average for WTTS ((Nh)w = 
1.8 x 10 21 cm" 2 ). 

The jVh values found from our spectral fits are roughly con- 
sistent with the visual extinctions Ay and the infrared extinc- 
tions Aj if we assume a standard gas-to dust ratio (Nn/Ay = 



Fig. 13. r av as a function of Nh for our sample. 



2 X 10 21 cm" z m ag' 1 , N H /A 3 



7.1 x 10 21 cirT 2 mag -1 ; 
Vuong et al.l |2003» . For a detailed discussion on the gas-to- 
dust ratio in TMC, we refer the reader to Glauser et al. (2007, 
in preparation). High photoelectric absorption might influence 
the spectral fits to low-resolution spectra, because the coolest 
plasma components, more affected by absorption, cannot be re- 
liably quantified. 

We studied possible biases introduced by high absorption 
by correlating Nh with T m and Lx- We have found Lx to range 
between approximately 10 28 erg s _1 and 10 31 erg s _1 , indepen- 
dent of the photoelectric absorption. The uncertainty of the de- 
termination of Lx, on the other hand, does increase with in- 
creasing Nh (and decreasing number of counts in the spectrum) 



as derived by Giide l et al. (2006a). Similarly, if we correlate T av 
with Nh (Fig.fTJll. we find a larger range of T m (symmetrically 
around log!T av » 7.0 [K]) for highly absorbed sources, while 
r av is found to be similar for all the sources with low absorp- 
tion (N H < 10 21 cirT 2 ), log T av * 7.0+0.2 [K]. The larger range 
of T av at higher Nh is likely to be the result of larger scatter due 
to less reliable spectral fitting. In any case, there is no trend 
toward higher r av for higher Nh- 

4.2.3. Correlation of L x with electron temperature 

In Fig. [14] we plot T av as a function of Lx and as a function 
of the X-ray surface flux (Fx) for CTTS and WTTS, respec- 
tively. The surface flux es have been calcul ated using the radii 
reported in Table 10 of lGiidel et a l. (2006a). In the plots for the 
WTTS we also show values for six m ain-sequence G-type so- 
lar analog stars (Telleschi et al., 2005), for 5 K-type stars (AB 
Dor from lSanz-Forcada et al.l2003l an d e Eri, 70 Oph A&B, 36 
Oph A&B from lwood & Linsky|20 06). and for 6 M-type main- 
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sequence stars (EQ Peg, AT Mic, AD Leo and EV Lac from 
Robrade & Schrmtjfcook AU Mic fro m lMagee et alJ2003l and 
Proxima Cen from lGiidel et alJl2004 . 

Spectral fits to low-resolution spectra that are subject to 
photoelectric absorption tend to ignore the coolest plasma com- 
ponents, as the soft part of the spectrum is most severely af- 
fected by the absorption. CTTS are on average more absorbed 
than WTTS. Given the larger absorptions, there could be a bias 
toward higher average temperatures in CTTS, although such a 
trend is not visible in Fig. [13] We nevertheless counteract a pos- 
sible residual bias by restricting the stellar sample used for the 
correlation to stars with Nh < 3 x 10 21 cirT 2 . The logarithmic 
means of jVh for CTTS and WTTS after these restrictions are 
1.2 X 10 21 cirT 2 and 8 x 10 20 cirT 2 , respectively, making these 
samples very similar with regard to absorption properties. 

Further, we exclude the very faint sources (with less than 
100 counts collectively in the three detectors) that could also 
produce unreliable Lx and L av results. In Fig.[l4]the filled cir- 
cles represent the stars used for the linear regression fit, while 
the stars excluded from the fit are plotted with small crosses. 
Overall, the absorbed and faint sources fit well to the trends 
found from less absorbed and more luminous sources, but their 
scatter tends to be larger. 

For CTTS, we find almost no correlation between L av and 
Lx or Fx (the correlation coefficients are 0.06 and 0.11, re- 
spectively). On the contrary, for WTTS L av is clearly corre- 
lated with both Lx and Fx. The correlation coefficients are 
0.69 and 0.72 for Lx and Fx, with 33 and 32 data points, 
respectively (Table Q~|i. The probability that no correlation is 
present is < 1% in either case. We computed the linear re- 
gression using the bisector OLS algorithm (no a priori re- 
lation between the two measured variables assumed) to find 
logL av = (0.23 + 0.03) log L x + (0.13 + 0.93) and logL av = 
(0.26 + 0.03) logLx + (5.21 + 0.24). WTTS follow a trend sim- 
ilar to that shown by MS stars in the L av vs. Lx relation. In the 
L av vs. Lx relation, on the other hand, we find that WTTS are 
in general hotter than MS stars for a given Fx- 

We have checked these results using the EM algorithm, 
finding slightly shallower slopes (see Table Q]). Shallower 
slopes are expected in th e EM algorithm wh en compared to the 
bisector OLS algorithm dlsobe et al.Ul990h . For CTTS, where 
no correlation is found, the two algorithms result in completely 
different sl opes (Table [D, an indication of absence of a linear 
regression dlsobe et aljl9 90; the different slopes in the absence 
of a correlation are a consequence of the defining minimization 
of the algorithm. The EM algorithm returns a slope of » 0, 
whereas the bisector algorithm yields a slope around unity). 

We use the Kaplan-Meier estimator in ASURV to compare 
the distributions of L av for CTTS and WTTS. The two distribu- 
tions are shown in Fig.[l5]for the case where we do not apply 
a restriction to A^: the solid line represents CTTS, while the 
dotted line represents WTTS. Only stars with more than 100 
counts in the three EPIC detectors are used. The probability 
that the two distributions arise from the same parent popula- 
tion is 0.7%-2%. If we restrict the sample to stars with low 
N H (< 3 x 10 21 cm" 2 ), we find the mean value of log L av [MK] 
= 7.10 [MK] with o- = 0.22 for CTTS and log L av [MK] = 6.88 
with cr — 0.17 for WTTS. The distribution is similar to the one 



found in Fig. [B] and the probability for the CTTS and WTTS 
distributions to originate from the same parent population is 
again only 0.7% -2%. 

We have checked if the difference found in the plasma tem- 
peratures of WTTS and CTTS could be attributed to abundance 
anomalies that may not have correc t ly been accounted f o r in th e 
spe ctral fits. iKastner et all d2002l) . IStelzer & Schmitll d2004 . 
and Drake et al J d2005h found large Ne/Fe and Ne/O abundance 
ratios in the spectrum of TW Hya. We have therefore fitted the 
spectra of the 19 CTTS with N H < 3 x 10 21 cm -2 and more than 
100 counts in the combined EPIC spectra (filled black bullets 
in Fig. [141 adopting abundances as found in TW Hya (O = 
0.2, Ne = 2.0, Fe _ 0.2, with respect to t he solar photospheric 
abundances of lAnders & Grevesse 1989, all other abundances 
as given in Sect. |3). The average temperatures obtained with 
this model are generally consistent within 0. 1 dex with the tem- 
peratures found based on our standard abundances. Only for 
one star, HO Tau AB, did we find an average temperature sig- 
nificantly lower, while the general trend toward higher L av for 
CTTS remains unchanged. We can therefore exclude that the 
difference in temperatures is induced by abundance anomalies 
as t hose observed in TW H ya. 

Telle schi et al. (2006a) derived the thermal structure of nine 
pre-main sequence stars from XEST based on high-resolution 
Reflection Grating Spectrometer data, using variable abun- 
dances. They found L av to be compatible with values used here, 
which were derived from EPIC CCD spectra (an exception is 
the CTTS SU Aur, for which the temperature found with RGS 
is even higher than that derived from the EPIC spectra). In the 
latter work, however, a difference in the abundances has been 
found between stars of spectral type K and stars of spectral type 
G. The abundances found for the K-type stars reflect approxi- 
mately the abundances used for the XEST EPIC fits (following 
an inverse FIP effect), while G-type stars show lower Ne/Fe 
and O/Fe abundance ratios. We therefore fitted the four G-type 
stars in our stellar samp le with an abundance p attern as found 
for this spectral class bv lTelleschi et al.ld2006al) . Again, we did 
not find a significant change in temperatures. In summary, we 
do not find any appreciable effect that abundance anomalies 
other than those adopted in our study might have on the tem- 
perature determination. We also note that Scelsi et al. (2006, 
in preparation) studied the abundances derived from the EPIC 
spectra of the brightest sources in the XEST sample and found 
average abundances very similar to the standard abundances 
used in our CCD fits. 

The difference in the coronal temperatures of CTTS and 
WTTS is in particular due to the larger L av found in the CTTS 
with low Lx- We therefore calculated the mean of all L av values 
for CTTS and WTTS with L x < 3 x 10 29 erg s 1 , low absorp- 
tion (A^h < 3 x 10 21 cirT 2 ), and more than 100 counts in the 
three detectors. For these stars, we computed the errors in L av 
as follows: We determined the 68% confidence contour on the 
y3 - To plane for these two "parameters of interest", i.e. the loci 
for which a fit can be achieved whose x 1 is larger by A;^ 2 = 2.3 
(1 cr) than the x 1 of the best fit- We then found the minimum 
and the maximum L av for this subset of models, and thus de- 
fined the error range for L av . Using these errors, we computed 
the weighted mean of log L av . For CTTS, we neglected DD Tau 
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Fig. 14. Left pannel: r av as a function of Lx for CTTS (top) and WTTS (bottom). Right panel: T av as a function of Fx for CTTS 
(top) and WTTS (bottom). The low-absorption TTS samples are marked by filled black bullets, while small cro sses give loci of 



high- absorption objects or sources with few counts (see text for details). Black diamonds mark solar analog stars (ITelleschi et al 



20051) and blue and green diamonds mark K- and M-type main-sequence stars, respectively (see text for references). The straight 
lines in the WTTS plots are linear regression fits (based on bisector regression, the dashed lines illustrating the error ranges in 
the slopes). 



AB, which shows extraordinarily high temperatures in two dif- 
ferent observations (two bullets at the hottest temperature in 
Fig.Qj). For CTTS, we find (log r av > = 6.97+0.06 (6.98+0.06 
if DD Tau AB is also considered) while for WTTS we find 
(log T av ) = 6.81 ± 0.05. We therefore conclude that the CTTS 
and WTTS with Lx < 3 x 10 29 erg s _1 are different at a 3<x level, 
fully supporting the significant differences in the regression fits 
that are based on the entire Lx range. 

In conclusion, we find the CTTS X-ray sources to be hotter 
than WTTS at a confidence level of > 98%, and this result is 
partly due to the presence of a Lx - T- dv relation for WTTS 
but its absence in CTTS. Further, the WTTS relation coincides 
with relations valid for main-sequence stars of different spectral 
types, including saturated and non-saturated stars at different 
evolutionary stages. 



5. Discussion 

5.1. Summary of trends 

We now discuss the trends and correlations described in the 
previous section and will also put them into a context with pre- 
vious reports, in particular from the COUP project. 




Fig. 15. Kaplan-Meier estimator for T av for CTTS (solid) and 
WTTS (dotted). Only sources with more than 100 counts in the 
three EPIC detectors are used. 

The most significant correlations that we reported above are 
those between stellar mass and Lx (slope ss 1 .7), between stel- 



14 



A. Telleschi et al.: X-ray emission from T Tau stars and the role of accretion 



lar bolometric luminosity L, and Lx (slope ~ 1), and between 
Lx and average electron temperature r av (slope 0.15-0.23), the 
latter applying only to WTTS. 

Further, we have found that Lx and Lx/L* are both lower, 
on average, for CTTS than for WTTS, each by a factor of » 2, 
compatible with the finding that the distributions of L, are sim- 
ilar for the two samples. In contrast, T av is, on average, higher 
by a factor of » 1 .7 for CTTS than for WTTS . 

Finally, we have studied possible correlations between 
Lx, T aw , or Lx/L* and the accretion rate but found at best un- 
convincing correlations. The same is true for a trend between 
L x and age. 

5.2. Comparison with previous studies 

The Lx - mass correlation. This relation has been reported 
prior to XMM-Newton and Chandra studies of star-forming 



region s, but with largely varying regressions. iFeigelson et al 



( 1993) found a slope of 3.6 + 0.6 for a sample of low-mass stars 
in the Chamaeleon I dark cloud based on ROSAT observations. 
There may be problems with more numerous upper limits at 
the low-mass (and low-l uminosity) end of the distribution in 
this study, as noticed by IPreibisch et al.l |2005). On the other 
hand, the COUP sample reveals a very similar correlation to 
ours, with a slope only marginally smaller (1.44+0.10) than for 
XEST (1 .69+0.1 1). The TMC sample thus essentially confirms 
the COUP results, and the residual difference might be due to 
a somewhat different distribution of stars in the HRD, perhaps 
indicating a different age distribution, as suggested from our 
discussion of this relation below. 

A clear difference between the two studie s is seen in the 
scatter around the regression curves. While Preibisch et al.l 



(2005) report a standard deviation of 0.65 dex (factor of 4.5) 
around the best-fit line, we find for the XEST sample values be- 
tween 0.38 (factor of 2.4, for WTTS) and 0.45 (factor of 2.8, for 
CTTS). We do not have a clear explanation for the smaller scat- 
ter in XEST, but note that i) similar findings apply to other cor- 
relations discussed below, and ii) the scatter found in the XEST 
results is close to the intrinsic uncertainties of any Lx measure- 
ment of magnetically active stars as these commonly va ry by 
such factors on various time scales: IPreibisch et al.l (12005) give 
a factor of 2 variation on long (yearly) time scales for the Orion 
sample. We are thus confident that the quality of our mass - 
Lx correlation corresponds to the minimum scatter that must 
be expected from snapshot observations of magnetically active 
stars. 

The Lx - L* correlation. The linear correlation between 
these variables expresses the classical result of X-ray satura- 
tion that has been found empirically for main-sequence stars 
( Vilhu & Rucinskii 1983 ). A similar law applies to v ery ac- 
tive m ain-sequence and subgiant stars (see review by iGtidell 
20041 and references therein), and certainly al so to pre-main se 



quence stars at various evolutionary stages (Fla ccomio et al 



Preibisch et al.l d2005l) find a well-defined linear correlation 
for WTTS (standard deviation around best-fit regression of 
0.52 dex), while for CTTS the scatter dominates (standard de- 
viation = 0.72 dex) and the relation is significantly flatter. The 
CTTS data points span a range of 3 orders of magnitude at 
a given L». In XEST, the standard deviation of the scattered 
points is only « 0.4 dex for CTTS, WTTS, and the entire sam- 
ple, with a range of Lx values at a given L„ of about 1 .5 dex. We 
are again not in a situation to explain the much tighter correla- 
tions for the XEST survey, but note that our spectral-fit method- 
ology may suppress numerical uncertainty introduced by pho- 
toelectric absorption that suppresses evidence of plasma com- 
ponents at lower temperatures. Because we used an emission- 
measure distribution with a prescribed low-temperature shape 
as usually found in magnetically active stars, the presence 
of the coolest components is interpreted based on the pres- 
ence of well-detected hotter plasma. An error analysis for Lx 
based in particular on A^h shows that for 76% of the sources 
in XEST, the intrinsic error range in Lx due to A^h is smaller 
than a factor of 3 (0.5 dex, in fact mostly much smaller), 
and the largest errors are o btained for faint sou rces subject to 
N H > (2 - 3) x 10 21 cm" 3 dGiidel et all l2006ah . We note that 
IPreibisch et al. ( 2005 ) used an X-ray luminosity averaged over 
the 10 days of exposure by Chandra, while in the XEST sam- 
ple we ha ve neglected time inter vals containing obvious flares. 
However. IPreibisch et al.l d2005l) found that the average Lx and 
the quiescent ("characteristic") Lx differ by a median factor of 
0.78, which would not be sufficient to explain the large scatter 
found in their Lx - L„ correlation for CTTS. 



2003b). Again, the COUP study is in complete agreement with 
our results, its regression slope being 1.04 + 0.06. There is, 
however, a significant difference between the XEST and the 
COUP studies once CTTS and WTTS are treated separately. 



Co mparison with main-sequence stars. IPreibisch et al 
(2005) also compare their Lx - L„ correlation with field main- 
sequence stars, and find a much shallower slope for the latter, 
but also a very large scatter. They similarly compare the Lx— 
mass relations with field stars, but there, they find a similar 
slope in the regression. We will not perform this comparison 
here, for the following reason. Field stars are found at various 
evolutionary stages, and as a consequence of stellar spin-down 
with age, the X-ray activity is subject to an evolutionary decay. 
Solar analogs decrease in Lx by three orders of magnitude from 
the z ero-age main sequence to the end of th e main-sequence 
life dGiidel et all 1 1 9971 iTelleschi et allEjoll) . Further, the evo- 
lutionary speed is different for G stars and low-mass M dwarfs, 
the latter remaining at relative ly high activi ty levels for a longer 
time (see Figs. 40 and 41 in lGiidelll2004h . Much of the scat- 
ter in Lx for a given mass or a given L„ is thus due to mass- 
dependent evolutionary decay, and any trend in Lx vs. L* de- 
pends on the stellar age distribution. In contrast, both in active 
main-sequence stars and TTS, no evolutionary effects are ex- 
pected for the Lx - L* relation if the stars are in a saturated 
regime, and therefore Lx °c L„. 

For the Lx - mass relation, the scatter in Lx for a given 
mass is only about one order of magnitude; this scatter is in- 
directly due to the scatter in L„ in the sample, due to different 
ages of stars of similar mass that contract vertically along the 
Hayashi track, provided that the X-ra y emission remain s in a 
saturated state (see HRD in Fig. 1 1 in lGiidel et al.ll200 6a). The 
scatter in Lx due to evolution on the main sequence is much 
larger (3 orders of magnitude) and is due to intrinsic decay of 
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the dynamo due to stellar spin down when the X-ray emission 
is no longer in a saturated state. 

The L x - T. dw correlation. A dependence between coro- 
nal electron temperature and emission measure or Lx (or nor- 
malized quantities such as the specific emissi on mea s ure or 
surfa ce X-ray flux ) has f irst been noted by Vaianal ( 1983 ) 
and ISchrijver et al. (11984. Quantitativ ely, for solar analogs, 
L x oc r 4 5±Q 3 dGiideletalJll997l see iGiidell Eool for a re- 
view ). For the pre-main s equence sample in the COUP sur- 
vey, iPjejbischetalJ (I2005D report a steep increase of the X-ray 
surface flux with the hotter temperature of their 2-component 
spectra, namely Fx <* T%, On the other hand, they find a rela- 
tively constant lower temperature, namely T\ m 10 MK. In our 
study, we apply a more physically appropriate continuous emis- 
sion measure distribution that does not distinguish between two 
isothermal components but that shows two power-law slopes 
on either side of the peak. The distribution of the logarithmi- 
cally averaged temperatures (Fig. [T3J does not show a pre- 
ferred value but a smooth distribution in the range 4-30 MK 
around a mean of 7.6 MK for WTTS and 12.6 MK for CTTS 
(Sect. 14.2.31 ). Our regression curve for the Lx - T- dv indicates 
Lx oc r 4 v 3 ~ 6 - 7 and Fx oc T dw & ~ 5 - 6 , compatible with the solar- 
analog relation as well as with the COUP relation for Ti. We 
note, however, that we find this trend only for WTTS (no sepa- 
rate analysis was provided for COUP). 

log{Lx/L*) distributions. Our distributions show the frac- 
tional X-ray luminosity of CTTS to be suppressed by a fac- 
tor of w 2 compared to WTTS, with a mean log(Lx/L*) ~ 
-3.39 + 0.06 and -3.73 + 0.05 for WTTS and CTTS, respec- 
tively. These values agree excellently with those of COUP: 
log(L x /L») ~ -3.31 and -3.7 3 for WTTS and CTTS, re- 
spectively dPreibischetal.ll2005l) . XEST contrasts with COUP 
in that CTTS are less X-ray efficie nt for all considered mass 
bins, whereas lPreibisch et al.l (120051) found no difference in the 
ranges 0.1-0.2 M Q and 0.5-1 M Q . 



5.3. The L x - mass relation 

Among the clearest correlations we have identified is the Lx- 
mass relation that closely corresponds to the finding in the 
COUP study. We now test the following: If we assume that 
TTS are in a saturated state (i.e. Lx oc L*, Fig. [5) and a re- 
lation between L» and stellar mass exists (related to the age 
distribution of the stars and details of the evolution along the 
pre-main-sequence tracks), then the relation between L x and 
mass could simply be a consequence of these two relations. 
Main-sequence stars follow the well-known mass-bolometric 
luminosity relation, which for stars in the mass range of 0.1- 
1.5M reads L* oc M 30 (from a regression analysis using the 
Siess et al. 2000 ZAMS data). For pre-main sequence stars, a 
mass-bolometric luminosity relation is not obvious; during the 
contraction phase, a star of a given mass decreases its L* by 
up to 2 orders of magnitude. However, if most stars in a sam- 
ple show similar ages, then an approximate mass-L» relation 
may apply to the respective isochrone. Fig. [16] illustrates the 
measured relation between L t and mass. The relation is rather 
tight, with a correlation coefficient of 0.90 for 113 sources. The 
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Fig. 16. Stellar bolometric luminosity as a function of mass. 
The straight line is a YfX OLS regression fit; dashed lines il- 
lustrate the 1-sigma errors in the slope. 



Y/X OLS regression gives log LJLq = (1.49 ± 0.07) log M + 
(0.23 + 0.04), with a standard deviation of 0.31. This relation 
can be compared with the theoretical prediction for an aver- 
age isochrone appropriate for the XEST sample. We found 
that the logarithmically averaged age of our targets is 2.4 Myr. 
Adopting the Siess et al. (2000) isochrones, we find, from a 
linear regression fit, a dependence L, oc M 124 , i.e., similar to 
the observed dependence and thus supporting our interpreta- 
tion. We note that the XEST sa mple is of course n ot located 
on an isochrone (see Fig. 11 in iGiidel et al]|2006al) . and that 
other evolutionary calculations may lead to somewhat different 
slopes of the isochrones. 

Adopting \og{Lxl L„) = -3.5 for our entire TTS sample 
(see Fig. [5}, we infer a relation logLx = 1.49 log M + 30.31, 
similar to the correlation found in Sect. 14. 1 . 1 1 In Fig. [17] we 
plot Lx as a function of mass after normalizing the observed 
Lx with Lx predicted from the above formula. The correla- 
tion found in Fig. Q] now disappears completely. The scatter 
in Fig.[l7]is due to the scatter in Fig. [5] and in Fig. [16] i.e., in- 
trinsic scatter not due to L„ or M, but for example due to the 
evolutionary decrease of L* and Lx along the Hayashi track for 
a given mass. 

We therefore conclude that the Lx-mass relation is not an 
intrinsic relation but a consequence of an approximate mass- 
luminosity relation for stars with similar ages, combined with 
a saturation law. 



Comparing th e stars in the HRD of the TMC (Fig. 11 in 
lGudelet alJl2006ah with t he HRD of the ONC sample (Fig. 
1 in iPreibisch et al. I l2005l) . we note that the ONC sample is 
somewhat younger, although it is not tightly arranged along 
one isochrone, but tends to show a somewhat flatter slope than 
the TMC sample. Th e slightly shallower Lx-mass relation in 
Preibisch et al. d2005l) may therefore be a straightforward con- 
sequence of the younger average age of the ONC sample. 
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Fig. 17. Lx as a function of mass after renormalizing Lx with 
the expected Lx based on mass-bolometric luminosity relation 
and the saturation law (see text for details). Different symbols 
mark different types of stars. 

5.4. Origin of the X-ray emission: CTTS vs WTTS 

The question that we address in this section is on the origin of 
the X-ray emission. How do CTTS and WTTS differ, and what 
may be the causes? Where and how are the X-rays formed? 

The relevant relations we have identified in this paper are 
the following: i) CTTS show, on average, a smaller X-ray lu- 
minosity in the EPIC band; ii) CTTS also reveal a significantly 
lower fractional X-ray luminosity, Lx/L*, than WTTS. iii) The 
average electron temperature in WTTS correlates with the to- 
tal L x for WTTS, but this is not the case for CTTS; the CTTS 
temperatures are on average significantly higher than those of 
WTTS. 

5.4.1 . Evidence for coronal emission 

The bulk of the X-ray emission described in this paper and in 
the COUP survey is consistent with an origin in a magnetic 
corona. The overall temperatures measured in TTS and the X- 
ray luminosities are similar to values also found in extremely 
active main -sequence and subgiant stars (e.g., review by 



Stel zer et all 12006c iFranciosini et all I2006I) clearly points to 



Gudelll2004l) Also frequent flaring in TTS dWolk et all I2005 



a coronal (or magnetospheric) origin of the X-ray emission. 
Although many of the T Tauri stars in the XEST sample are 
thought to be fully convective and hence unable to support a 
solar-like <?-Q dynamo, fully-convective main-sequence stars 
do show magnetic activity, and dynamo mec hanisms have been 
proposed that may oper ate in such stars (e.g.. | Dobler et al .2006 
and references therein; Kii ker & Riidigei 1999L . A discussion 
of what XEST can tell us about t he dynamos acting i n T Tauri 
stars in Taurus-Auriga is given by Bri ggs et al.l (|2006). 

We have found in the previous section that the Lx - mass 
relation is due to saturation and a mass-bolometric luminosity 
relation for TMC stars. The most fundamental relation is in- 



deed the one between Lx and L„, which indicates that the PMS 
stars in Taurus are saturated at log(Lx/L») a -3.5, in analogy 
to main-sequence stars. This result suggests that the bulk of X- 
ray emission in PMS stars by CCD detectors used in XEST and 
in COUP arises from a magnetic corona, and we therefore sug- 
gest that similar processes as in MS stars should be responsible 
for the dominant energy output from the hot plasma (we will 
rediscuss emission from the softest X-ray emitting plasmas be- 
low). 



5.4.2. Is a lower L x intrinsic to CTTS coronae? 



We have tested (Sect. I4.2.U whether the lower Lx seen in 
CTTS compared to WTTS could be a consequence of gener- 
ally lower L, of CTTS, provided that all stars are in a similar 
saturation regime. This can be explicitly rejected because we 
found the L* distributions for the two samples to be similar, 
with a trend that it is rather the CTTS sample that is slightly 
more luminous. Also, the Lx/L, distributions are explicitly dif- 
ferent, offset by about the same factor as the Lx distributions 
themselves. 

In previous works the lower X-ray activity of CTTS com- 
pared to WTTS in Taurus-Auriga has been attributed to the 
slower rotation of CTTS and an anticorrelation of activity and 
rotation period as exh ibited by active solar-like stars (e.g., 
iNeuhauser et al.l Il995b . However, the lower X-ray activity of 
CTTS has been observed in other star-forming regions where 
an anticor relation of activity an d rotation period is cl early not 
seen (e.g. IPreibisch et alJl2005h . iBriggs et al.l |2006) demon- 
strate that an apparent activity-rotation relation in Taurus- 
Auriga naturally results from the dependences of activity on 
mass and accretion status reported here and in other star- 
forming regions because the fast rotators in Taurus-Auriga are 
mainly higher-mass and non-accreting while the slow rotators 
are mainly lower-mass and accreting. There is no convincing 
evidence for an anticorrelation of X-ray activity and rotation 
period in T Tauri stars, and therefore no evidence that the lower 
activity of CTTS is due to their slower rotation. 

Indeed, even if a solar-like dynamo operates in T Tauri 
stars, their long convective turnover timescales lead to the 
expectation that all stars with measured rotation periods 
s hould have saturat ed (or supersaturated) emission (e.g., 



Preibisch et al.ll2005l) and show no anticorrelation of activity 
and rotation period. 

Different internal structure in WTTS and CTTS (unless in- 
duced by accretion processes) is not a likely explanation either. 
As was recognized in early infrared and optical surveys of the 
TMC stellar popul ation, WTTS and CTTS occ upy the same 
region in the HRD dKenvon & Hartmannl ll995), with no evo- 
lutionary separation, indicating that the transition from CTTS 
to WTTS occurs at very different ages for different stars. 

This suggests that the lower Lx of CTTS is an intrinsic 
property of a corona that is heated in the presence of an ac- 
cretion disk and active accretion onto the star. We now briefly 
consider the influence of accretion on coronal heating. 
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5.4.3. The role of active accretion 

Accretion has variously been suggested to enhance or to 
suppress plasma heating. First, accretion hot spots may 
heat plasma to temperatures in excess of one million K 
as gas shocks near the surface in nearly free fall (e.g., 
Calvet & Gullbringll998 ). High-resolution X-ray spectroscopy 
has provided some indirect evidence that such accretion- 
induced X-rays mi ght constitu t e an im portant part of the mea- 
surable spectra. iKastner et al.1 (12002b have interpreted excep- 
tionally high-densities and very cool (T ~ 3 MK) X-ray 
emitting material in the CTTS TW Hya as being the re- 
sult of accretion shocks. This idea was further elaborated by 
1 J k 

Stelzer & Schmitt ( 2004) who also suggested that anomalously 
high Ne and N abundance in this X-ray source indicates that re- 
fractory elements such as Fe condense onto dust grains in the 
disks, and that this material is eventually not accreted. A small 
number of additional X-ray spectra have been studied, but the 
situa tion is complex and con tradictory: Schmitt et al.1 (2005) 
and Robrade & Schmitt (2006) report intermedi ately high den- 



sities in BP Tau (« e ~ 3 x 10 cm ), and Giint her et al 



(2006) 1 ind similarly high densities in V4046 Sgr. On the other 
hand, very low densities an d no strong abundan ce anomalies 
have been found in T Tau ( Giidel et al. . 2006cl) a nd the ac- 
creting Herbig star AB Aur (Tell eschi et all l2006bh . However, 
Telleschi et al.l (I2006al) found evidence that CTTS in general 
maintain an excess of cool (1-4 MK) plasma compared to 
WTTS. This expresses itself in Ovu line fluxes that are sim- 
ilar to the flux in the O vm Ly<? line, a condition very different 
from WTTS where the O vn lines are very faint (similar to mag- 
netically active ZAMS stars). 

This soft excess seems - as far as the still small statis- 
tical sample of stars suggests - to be related to accretion in 
CTTS. But where are the X-rays produced? For T Tau and AB 



Aur, a production in accretion shocks is unlikely ( Gii del et al 
2006c; iTelleschi et all l2006bh while accretion shocks may 



be responsible for the softest emiss i on in TW Hya, BP Tau, 
and V4046 Sgr dKastner et all 120021: IStelzer & Schmitl 12004 



Schmitt et all 12005; Giint her et all 2006). 



5.4.4. Coronal modification by accretion streams? 

An alternative possibility is that accretion influences the mag- 
netic field structure of the stars, or that the accreting material is 
changing the h eating behavior in the coronal magnetic fields. 
Preibisch et al.l (120051) proposed that the mass-loaded loops of 
accreting stars are denser than the loops of WTTS, so that when 
a magnetic reconnection event occurs, the plasma would be 
heated to much lower temperatures outside the X-ray detection 
limit. This could account for the deficiency of X-ray luminos- 
ity in CTTS compared to WTTS . A similar scenario has been 
proposed by lGiidel et al.l (I2006cl) to specifically explain the soft 
excess observed in T Tau. In this scenario, the magnetospheric 
geometry is influenced by the accretion stream. A fraction of 
the cool accreting material enters the coronal active regions and 
cools the magnetic loops there for three reasons: i) the accre- 
tion flow may reorganize magnetic fields, stretching them out 
and making them less susceptible to magnetic reconnection; ii) 



the accreting material is cold, lowering the resultant tempera- 
ture on the loops when mixing with the hot plasma; and iii) the 
accreting material itself adds to the coronal density, inducing 
larger radiative losses and more rapid cooling. These result- 
ing soft X-rays from cool material can be detected by the RGS 
instruments that are sensitive at low energies and provide the 
spectral resolution to record flux from lines formed at tempera- 
tures below 3 MK, such as O vn and N vi, but the same emission 
will not be separately identified by the EPIC detectors that are 
relatively insensitive at the rele vant temperature, an d provide 
only very low energy resolution. Giide l et al. (2006c) estimated 
for T Tau that only of order 1 % of the accreting material would 
be needed to penetrate active regions on the star and be heated 
to 2 MK. 

If the total coronal energy release rate (averaged over time 
scales longer than energy release events such as flares) is de- 
termined by the stellar dynamo that forms a magnetic corona 
(and by convective properties near the stellar surface), then we 
would expect similar coronal radiative losses for CTTS and 
WTTS. Could it be that the radiative output from the corona 
is indeed equivalent but that part of the coronal emission has 
shifted to the softest part of the spectrum, remaining undetected 
in EPIC CCD spectroscopy while detected as a soft excess by 
RGS? At least in the case of T Tau, the 0.3-10keV X-ray lumi- 
nosity has been severely underestimated by EPIC CCD analysis 
alone, as these instruments missed the softest component, also 
subject to considerable photoelectric absorption, and suggested 
Lx to be only « 6 0% of Lx determine d from the combined 
RGS+EPIC spectra dGiidel et alll2006ch . 
We systematicall y studied our RGS 



spectra (from 

Telleschi et alJ 2006ah to find out whether the soft excess in 



our CTTS sample provides the "missing luminosity". We have 
therefore compared the unabsorbed 0.1-0.5 keV X-ray lumi- 
nosity from the EPIC spectral fits ( Giidel et all l2006al) with 
Lx in the same range from the combined EPIC+RGS fits 
( Telleschfetall |2006a). The comparison is useful only for tar- 
gets that do not suffer from strong absorption; this is the case 
for the CTTS BP Tau and DN Tau and the WTTS HD 283572, 
V773 Tau, and V410 Tau. We found no systematic difference 
to explain the factor of 2 underluminosity of CTTS. It appears 
that the EPICs record the soft emission from the coolest plasma 
sufficiently well to register similar Lx as the RGS detectors, 
but the temperature discrimination is clearly inferior to high- 
resolution spectroscopy. Also, the softest range is still domi- 
nated by continuum emission from hotter plasma, and the soft 
excess in these stars provides relatively little spectral flux. 

We have found, on the other hand, that CTTS show, on av- 
erage, higher electron temperatures (averaged over the compo- 
nents detected by the EPIC cameras) than WTTS. This could 
be an effect of depletion of the intermediate and cooler tem- 
perature ranges by the accretion process, as suggested above, 
thus moving the average temperature of the detected coronal 
components to higher temperatures. We have therefore tested 
whether the harder portion of the EPIC spectra which is ra- 
diated by the hottest coronal plasma components also shows 
a statistical difference between CTTS and WTTS. We chose 
the 1.5-10 keV range for this test. We found, however, the 
same discrepancy between the two stellar groups, suggesting 
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that the emission measures of the hottest plasma components 
themselves are also suppressed in CTTS compared to WTTS. 

We extended our comparison to the XEST results obtained 
from EPIC only, but again found CTTS to be underluminous 
by similar factors in different X-ray energy bands. To explain 
the deficit of X-ray emission in CTTS, it could thus be that the 
accretion process is cooling active-region plasma to an extent 
that it is also no longer detected in the RGS band. 

5.4.5. Coronal heating in WTTS and CTTS due to 
flares? 

This brings us to the correlation between average coronal tem- 
perature and Lx which is present in WTTS but absent in CTTS. 
WTTS show a trend in which T m increases with Lx, and this 
trend is the same as previously found for main-seq uence so- 
lar analogs dGiidel et all 1 1 997t iTelleschi et all 120051) . We note 



that this relation remains valid even for stars with different X- 
ray saturation limits and different L„ (see Fig. FBI. M dwarfs at 
the saturation limit reveal lower coronal temperatures than K 
or G dwarfs at their satur ation limit. The cause for this relation 
is not clear. iGiidell (120041) pointed out that the slope of the re- 
gression function (using emission measure instead of Lx) is the 
same as the slope of the regression between peak t e mpera ture 
and peak emission measure in stellar flares. Giidel (2004) hy- 
pothesized that coronal emission is formed by a superposition 
of continuously occurring "stochastic flares", with the conse- 
quence that larger, hotter flares that occur more frequently in 
more active stars not only produce the dominant portion of the 
observed emission measure, but also heat the observed plasma 
to higher temperatures than in lower-activity stars. The larger 
rate of large flares would be a consequence of denser packing 
of magnetic fields in more active stars, inducing more frequent 
explosive magnetic r econnection, includ ing larger flares than 
in low-activity stars ( Giidel et all 1997 ). A similar trend for 
WTTS as for main-sequence stars is therefore perhaps not sur- 
prising: X-ray production in both types of stars is thought to 
be entirely based on the magnetic field production by the in- 
ternal dynamo. This analogy fully supports solar-like coronal 
processes in WTTS. 

The Lx-Tw correlation is absent in CTTS. The distinguish- 
ing property of CTTS is active accretion, which thus is most 
likely the determining factor for the predominantly hot coronal 
plasma. Temperatures like those determined as r av in the XEST 
survey cannot be produced in accretion shocks, again pointing 
to a predominantly coronal origin of the hot, dominant plasma 
component. If the flare-heating concept has merit in CTTS as 
well, then it seems that flares in CTTS are predominantly hot, 
even if Lx is low. We can only speculate about the origin of 
this feature. A possibility are star-disk magnetic fields, but it 
is unclear why flares occuring in such loop systems should be 
hotter. 

6. Summary and conclusions 

We have studied X-ray parameters of a large sample of X-ray 
spectra of CTTS and WTTS in the Taurus Molecular Cloud. 
Our principal interest has been in a characterization of X-rays 



in the two types of stars, in finding correlations between X- 
ray parameters and fundamental stellar properties and among 
themselves, and most importantly in comparing our findings 
between accretors and non-accretors. This study has been mo- 
tivated by numerous previous reports on correlations and dif- 
ferences between CTTS and WTTS in nearby star-forming re- 
gions, and in particular by the COUP study of the Orion Nebula 
Cluster. The XEST project has provided the deepest and, for 
the surveyed area, most complete X-ray sample in the Taurus 
region to date. We have used a CTTS and a WTTS sample of 
comparable size. 

We have correlated Lx and the average coronal tempera- 
ture, r av , with various stellar parameters, and conclude the fol- 
lowing from our study: 

- The X-ray luminosity is well correlated with the stellar 
mass, with a dependence Lx M lJ , similar to what has 
been shown in COUP and previous TTS studies, but we 
find that this correlation is only an expression of satura- 
tion and a mass-(bolometric) luminosity relation for our 
pre-main sequence sample. As long as the stellar sample is 
saturated, Lx is a function of L„, and the latter is correlated 
with stellar mass for a given isochrone. From stellar evolu- 
tion calculations (e.g., Siess et al. 2000l) . the functionality 
between mass and L* can be derived. This is fully analo- 
gous to main-sequence stars where approximately L„ oc M 3 
holds. For a typical isochrone of TMC stars with ages of 
2-3 Myr, the exponent is smaller as can be seen on a pre- 
main sequence HRD (Fig. 11 in iGiidel et al.ll2006ah . For 
our sample, L„ oc M L5 . 

- A saturation relation holds for both CTTS and WTTS, al- 
though Lx/L* is, on average, smaller by a factor of 2 for 
CTTS compared to WTTS. 

- We find that the distributions of L„ are similar for CTTS and 
WTTS. As a consequence, we find a significant difference 
in the X-ray luminosity functions for CTTS and WTTS, 
the former being fainter by about a factor of two. The sup- 
pressed X-ray production in CTTS is thus intrinsic to the 
source and not due to selection bias. 

- We emphasize that the lower X-ray production in CTTS 
refers to the range of plasma temperatures accessible by 
CCD cameras such as those used here and in COUP. It is 
possible that some of the energy release is shifted to lower 
temperatures outside the range easily accessible to CCD de- 
tectors. Those soft regions of the X-ray spectrum are also 
subject to increased photoelectric absorption, which makes 
detection of cool plasma more difficult. CTTS are indeed 
more absorbed than WTTS, namely by a factor of « 2.5. 

- We investigated whether all X-ray spectral ranges show X- 
ray suppression in CTTS. The hardest portion (1.5-10keV) 
shows the X-ray deficiency in CTTS vs. WTTS indepen- 
dently, even though CTTS reveal higher average tempera- 
tures. We hypothesized that a fraction of the emission mea- 
sure has been cooled to poorly detectable or undetectable 
temperatures in CTTS. CTTS indeed show a soft excess 
in their high-resolution X-ray spectra, char acterized by un 



usual ly strong O vn lines from cool plasma (Tell eschi et al 



2006a). These lines cannot be resolved by EPIC. We there- 
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fore checked whether R GS sp e ctrosco py of little absorbed 
CTTS stars in iTelleschi et all d2006al) (DN Tau, BP Tau) 
indicates a relative increase of Lx in the soft 0.1-0.5 keV 
band relative to WTTS, but found no significant effect. The 
soft flux may have been sufficiently well detected by EPIC 
in these low-absorption stars (but with little temperature 
discrimination). Also, the softest range is still dominated 
by continuum emission from hotter plasma, and the soft 
excess in these stars provides relatively little spectral flux. 
The X-ray deficiency in CTTS thus remains. 



The situation is clearly different in T Tau ( Gudel et al 



2006ch : in this much more strongly absorbed source, a 



very large amount of very cool X-ray emitting plasma was 
detected based exclusively on anomalously strong Ovu 
line emission in the grating spectrum but went unnoticed 
in CCD spectroscopy. The analysis of the latter spectrum 
alone led to an underestimate of the 0.3-10 keV luminosity 
by 40%. 

In conclusion, it seems that the entire X-ray range acces- 
sible to CCD spectroscopy reveals suppressed X-ray emis- 
sion compared to WTTS, although additional components 
may be present at cool temperatures that may be missed by 
the CCD spectra, especially if A^h is sufficiently high. 

- A possible cause for the suppression of X-ray emission in 
CTTS may be the accretion streams themselves. If only a 
small portion of the accreting matter penetrates into hot 
coronal magnetic structures, the plasma may cool as more 
matter needs to be heated and as the increase in density 
incre ases the cooling effi ciency. This may lead to a soft ex- 
cess ( Giidel et al. , 2006ch . or to a co oling of plasma totem- 
peratures outside the X-ray regime (iPreibisch et all 12005). 
so that a significant deficiency of X-ray emission may be 
measured in the spectral range that is accessible to CCD 
cameras, and that is not subject to significant photoelectric 
absorption. The soft excess and the hot-plasma deficiency 
seem to be related to the presence of accretion. 

- X-ray production in shocks at the base of accretion streams 
has been suggested previously from high-resolution spec- 
tra of CTTS. The shocked plasma would add soft emis- 
sion to the spectra as well, but again, CCD spectroscopy 
may miss this emission, or the latter may be subject to 
absorption. Our CCD survey does not provide the appro- 
priate means to test X-ray production in accretion shocks 
in CTTS, and can therefore also not exclude such mecha- 
nisms. High-resolution grating spectroscopy is required. 
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